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ABSTRACT 

Using direct simulations of hydromagnetic turbulence driven by random polarized waves it is shown that 
dynamo action is possible over a wide range of magnetic Prandtl numbers from 10"^ to 1. Triply periodic 
boundary conditions are being used. In the final saturated state the resulting magnetic field has a large-scale 
component of Beltrami type. For the kinematic phase, growth rates have been determined for magnetic Prandtl 
numbers between 0.01 and 1, but only the case with the smallest magnetic Prandtl number shows large-scale 
magnetic fields. It is less organized than in the nonlinear stage. For small magnetic Prandtl numbers the growth 
rates are comparable to those calculated from an alpha squared mean-field dynamo. In the linear regime the 
magnetic helicity spectrum has a short inertial range compatible with a -5/3 power law, while in the nonlinear 
regime it is the current helicity whose spectrum may be compatible with such a law. In the saturated case, the 
spectral magnetic energy in the inertial range is in slight excess over the spectral kinetic energy, although for 
small magnetic Prandtl numbers the magnetic energy spectrum reaches its resistive cut off wavenumber more 
quickly. The viscous energy dissipation declines with the square root of the magnetic Prandtl number, which 
implies that most of the energy is dissipated via Joule heat. 
Subject headings: MHD - turbulence 



1. INTRODUCTION 

Many astrophysical plasmas are turbulent and tend to be 
magnetized. The magnetic fields can have typical length 
scales that are either larger or smaller than that of the energy- 
carrying eddies. We speak then correspondingly of large- 
scale or small-scale dynamos. Small-scale dynamos can al- 
ready work in statistically mirror-symmetric isotropic homo- 
geneous turbulence, whereas large-scale dynamos require in 
general a departure from parity-invariant or mirror-symmetric 
flows. The excitation conditions of small-scale dynamos de- 
pend sensitively on the value of the magnetic Prandtl number, 
i.e. the ratio of kinematic viscosity to magnetic diffusivity, 
Pim = v/77, where v is the kinematic viscosity and 77 the mag- 
netic diffusivity. This sensitivity is related to the fact that in 
the kinematic regime the spectral magnetic energy is peaked 
at the resistive scale. As was pointed out originally by Ro- 
gachevskii & Kleeorin (1997), and more recently by Boldyrev 
& Cattaneo (2004), the slope of the kinetic energy spectrum 
is important for the onset of small-scale dynamo action. It 
matters therefore whether the resistive scale lies in the viscous 
range (Pr^ ~ 1), within the inertial range (Pr^ 5 0.1), or right 
within the range where the bottleneck occurs (Pr^ ~ 0.1). 
The bottleneck effect refers to the spectral subrange just be- 
fore the dissipation range where the kinetic energy spectrum 
is shallower than in the inertial range. Within the bottleneck 
range the velocity increments diverge even more strongly with 
decreasing separation than in the inertial range, making dy- 
namo action harder still. Indeed, for small values of Pr^ 
the critical value of the magnetic Reynolds number above 
which small-scale dynamo action occurs increases therefore 
sharply toward Pr^r = 0.1 (Schekochihin et al. 2005), and 
then decreases slightly for Pr^ < 0.05 (Iskakov et al. 2007). 
However, even with the computing power available today, di- 
rect simulations of small-scale dynamo action are still only 
marginally possible at such small values of Pr^. 

For certain types of flows dynamo action is easier to achieve 
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even though Pr^ is small. The Taylor-Green flow is an exam- 
ple where the critical value of the magnetic Reynolds number 
becomes constant for Pr^ < 0.1 (Ponty et al. 2004, 2005). In 
this flow there can be large-scale patches with finite kinetic 
helicity of opposite sign. A completely different example is 
fully helical turbulence where the excitation condition for dy- 
namo action is virtually unchanged as Pr^ decreases from 
1 to 0.1 (Brandenburg 2001, hereafter BOl). For an ABC- 
flow dynamo, Mininni (2007) found a weak dependence of 
the threshold value of the magnetic Reynolds number Re^ on 
Pr^. For PrM < 0.1 the threshold value seemed to become 
asymptotically independent of Pr^ and dynamo action was 
demonstrated for values of Pr^ down to 5 x lO^^*. 

In many astrophysical bodies, Pr^ is indeed rather small 
(around 10"^). Such systems still possess dynamo action and 
can have large-scale magnetic fields. It is likely that such sys- 
tems belong to the second class of systems where the excita- 
tion conditions are not drastically altered toward small values 
of Pr^. Indeed, large-scale magnetic fields are found regard- 
less of whether Pr^ is small (e.g., the Sun and other stars 
with outer convection zones, as well as planets) or large (e.g., 
in spiral galaxies, because of their low densities). 

The purpose of this paper is to point out that a strong Pr^ 
dependence does not occur in systems where the magnetic 
field generation is predominantly due to a large-scale dynamo. 
Such systems have been studied in idealized settings such 
as periodic boxes using explicit forcing functions for driv- 
ing the turbulence. This has significant advantages in that 
periodic boundary conditions can be used, energy spectra are 
easily computed and, most importantly, isotropy and homo- 
geneity eases comparison with turbulence theory. A disadvan- 
tage is that the magnetic helicity can only change on resistive 
timescales, which slows down the saturation (BOl). 

With these provisions in mind, we consider now simula- 
tions of maximally helical turbulence in triply periodic boxes 
where we keep in most cases the fluid Reynolds number. 
Re - UnTis/vkf, constant and vary the magnetic Reynolds num- 
ber, RtM = u^rnf,Mkf, and thereby Pr^r (= Re^/Re). Here, 
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Mims is the rms velocity of the turbulence and kf is the forcing 
wavenumber. According to BOl the dynamo should be ex- 
cited whenever the domain is large enough (2-3 times larger 
than the forcing scale) and the magnetic Reynolds number 
exceeds unity (Re^ > 1.1... 1.4 or so). This was confirmed 
for magnetic Prandtl numbers as low as 0.1. In the present 
work we consider kinematic dynamo action down to values 
of Pim - 10"^ and nonlinear saturated dynamos down to 
PrM = 10-^ 

2. THE METHOD 

We solve the hydromagnetic equations for velocity U, log- 
arithmic density In p, and magnetic vector potential A for an 
isothermal gas in the presence of an externally imposed heli- 
cal forcing function /, 

ot 
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Here, B - V x A is the magnetic field, / = V x B/yUo is 
the current density, fio is the vacuum permeability, c, is the 
isothermal speed of sound, and S,j - ^(Uij+Ujj)--^Sij'V-U is 
the traceless rate of strain tensor We consider a triply periodic 
domain of size L^, so the smallest wavenumber in the domain 
is ki = 2jt/L. The forcing function consists of eigenfunctions 
of the curl operator with positive eigenvalues and is therefore 
fully heHcal with / ■ V x / - kf^, where 3.5 < k/ki < 4.5 is 
wavenumber interval of the forcing function, whose average 
value is referred to as ^/ - 4ki. The amplitude of / is such 
that the Mach number is Mims/cj ~ 0.1, so compressive effects 
are negligible (Dobler et al. 2003). 

The initial conditions consist of a weak Beltrami field. The 
initial velocity is zero and the initial density is uniform with 
p - Pq - const. Note that the volume-averaged density re- 
mains constant, i.e., (p) = po. 

The model is equivalent to that of BOl, except that there 
the value of kfjki was chosen to be either 5 or 30. In order 
for large-scale dynamo action to be possible, kj lki must at 
least be larger than 2 (Haugen et al. 2004), but 3 is already 
sufficient (Brandenburg et al. 2008). In the kinematic regime 
the fastest growing mode is expected to have the wavenumber 
kfjl, so in order that this wavenumber is distinct from ki, we 
have chosen kfjki — A throughout this paper. 

3. RESULTS 

We begin by presenting results for Re s; 670 where we vary 
PrM in the range 0.01 < Pr^ < 1, i.e., Re^ is varied in the 
range 6.7 < Re^ < 670. The dynamo is excited in all those 
cases, but the growth rate A varies. We consider first the kine- 
matic regime where the magnetic field is weak and turn then 
to the nonlinear regime where the magnetic field has satu- 
rated. Most of the results presented below have been obtained 
at a resolution of 512^ meshpoints. The solution was first 
evolved at lower resolution (128^ meshpoints), then remeshed 
to twice the resolution, again evolved for some time, and fi- 
nally remeshed to 512'' meshpoints, and again evolved for 
some time. Data for the kinematic regime are only used after 
the initial transients have disappeared and a clear exponen- 
tial growth has developed at all length scales for at least some 



40 turnover times (also for the runs with a resolution of 512^ 
meshpoints). The run with 128"^ meshpoints has been evolved 
all the way into saturation, and it was then remeshed twice by 
a factor of 2, just like in the kinematic regime. 

3.1. Field structure in the kinematic regime 

Visualizations of one component of the magnetic field show 
the emergence of a large-scale magnetic field for small val- 
ues of Pr^f. This is clearly demonstrated in Fig.[Tl where we 
see for Pr^ = 0.01 a large-scale pattern with a systematic 
variation in the y direction. For Pr^ =0.1 there is no such 
variation, although there are some extended patches in which 
the field orientation is the same. For Pr^ - 1 even this is no 
longer the case and the field appears completely random with 
small-scale variations only. 

We emphasize that random and patch-like structures only 
occur in the kinematic regime. In the saturated regime a large- 
scale field emerges in all cases. This will be discussed in ^3.71 

3.2. Growth rates 

The growth rate is calculated as the average of the instanta- 
neous growth rate, d In Bims/df. Examples are shown in Fig.|2] 
for runs with Re = 670 and different values of Pr^ using 512^ 
meshpoints. In Fig. [3] we show growth rates normalized by 
Ui-mskf (inverse turnover times) as a function of Re^ for three 
values of Rcjvf and compare with the corresponding results for 
non-helical turbulence forced at larger scales in the wavenum- 
ber interval 1 < k/ki < 2 (Haugen et al. 2004). In that case 
we use kf = 1.5 for the average value. The small-scale dy- 
namo is then only excited when Re^ > 35. For Re^ > 100 
the growth rates for helical turbulence with kf = 4 are quite 
similar to those of non-helical turbulence with kf - 1.5. We 
have also calculated growth rates for the non-helical case with 
kf = 4 and find the same values as in the helical case. We note 
that in all cases, and even for small values of Pr^, the growth 
rates based on the rms magnetic field are equal to those based 
on the rms values of the mean fields obtained by averaging 
over any two coordinate directions. 

In both helical and non-helical cases, when Re^ is large 

1 /2 

enough, A increases like Re^ , as expected (Schekochihin et 
al. 2004). This is because the eddy turnover rate at the re- 



sistive scale is oc k^^^, but because k,j/kf cc Re^"* we have 

1 /2 

A oc Re^ . However, when there is also large-scale dynamo 
action, one expects there to be a lower bound for A given by 
the growth rate for the large-scale dynamo, /Ils. Using the 
theory for an dynamo (Moffatt 1978, Krause & Radler 
1980), we have 
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^LS = \ak\ -{?] + r],)k^ 
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where a is a pseudo scalar (the a effect) and rjt is the turbu- 
lent magnetic diffusivity. For fully helical turbulence we have 
lal X Mi-ms/3 and rj, x u^-n^.l'ikf (Sur et al. 2008), so we can 
write the growth rate as 
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where we have put \k\ - k\ . Over the parameter range con- 
sidered in this paper (Re^ > 6.7 and kflk\ - 4), /tLs/wrms^/ 
increases only slightly from 0.053 to 0.062 as Re^ increases. 

The result shown in Fig. [3] gives values that are systemat- 
ically below /Ils. There could be two reasons for this dis- 
crepancy. On the one hand, the accuracy of the estimates 
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Fig. 2. — Instantaneous growth rate for runs with different values of Pr^ 
and 512' meshpoints. The straight lines give the growth rates obtained by 
averaging over the indicated time interval. The solid, dotted, and dashed 
lines are for Pr^ = 0.01, 0.1, and 1, respectively. 
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Fig. 3. — Dynamo growth rates of the rms magnetic field for helical turbu- 
lence with Re = 670 (solid line) compared with growth rates for non-helical 
turbulence (dashed lines; adapted from Haugen et al. 2004). 

~ Mrms/3 and 77, a: u^ms/^kf may not be good enough. On 
the other hand, equation ^ is only an approximation in cases 
where /Ils ^ 0, because then memory effects become impor- 
tant. This means that, when allowing a and rj, to be integral 
kernels in time, they are no longer proportional to 6 functions, 
but have finite widths in time. This effect has recently been 
studied by Hubbard & Brandenburg (2008) and can be quite 
dramatic in some cases. 

In Fig. [3] we have varied Rcm by changing Ptm and keep- 
ing Re - const = 670. We can therefore also consider this 
graph as a representation of the magnetic Prandtl number de- 
pendence of A. However, Pr^ can also be changed while 



Fig. 4. — Dependence of dynamo growth rates of the rms magnetic field on 
Pi'M for helical turbulence with Re^ = 6.7 (dashed line) and Re = 670 (solid 
line). Here the solid line corresponds to the solid line in Fig. [3] 

keeping Re^ = const = 6.7. The corresponding result is 
shown in Fig. |4] and compared with the previous case. The 
two graphs are in reasonable agreement for small values of 
Pr^, but the rise of A for Re = 670 around Pr^ = 1 is not 
seen in the case with Re^ = 6.7. This suggests that the tran- 
sition from a purely large-scale turbulent dynamo to a mixed 
large-scale and small-scale turbulent dynamo requires values 
of ReM above some critical value (somewhere between 10 and 
100), and is not just determined by the value of Pr^- 

3.3. Spectra 

The transition from a purely large-scale turbulent dynamo 
to a mixed large-scale and small-scale turbulent dynamo is 
accompanied by characteristic changes in the spectral proper- 
ties of the magnetic field. In the following we employ shell- 
integrated spectra of kinetic and magnetic energy, E(k) and 
M(k), respectively, as well as of kinetic and magnetic helic- 
ities, F{k) and H(k), respectively. These spectra are normal- 
ized such that / E{k)dk = \{U^} = E, j M{k)dk = \_{B^} = 

M, J F(k)dk = {W ■ U), and ] H{k)Ak = <A ■ B), where 
W = V X t/ is the vorticity. 

In Fig. |5] we plot E{k) and M{k) for three cases with 
PrM - 1, 0.1, and 0.01, keeping Re = 670 in all cases. The 
magnetic energy spectra are compensated by exp(-/lf), where 
A is the numerically determined growth rate for each run, and 
then averaged in time. We consider here only the kinematic 
regime when the magnetic energy is weak. The kinetic energy 
spectra are then always the same. Since the magnetic energy 
is weak, we have scaled the magnetic energy spectra for differ- 
ent PrM to an arbitrarily chosen reference value of 10"^ below 
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Fig. 5. — Spectra of kinetic and magnetic energies in the kinematic regime 
forPr^ = 0.01, 0.1, and 1. 
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Fig. 6. — Compensated spectra of kinetic and magnetic energies and helic- 
ities in the kinematic regime for Pi'm = 1 . The spectra are denoted by letters 
E, F, M, and H. as described in the text. 

the kinetic energy spectrum. For Pr^ = 1 the magnetic energy 
seems to follows an approximate Kazantsev (1968) spectrum 
with a range proportional to k^^^ and is peaked at the resistive 
scale near k/ki = 50. For smaller values of Pr^ the peak of 
magnetic energy moves to smaller wavenumbers. 

In order to judge the correspondence with various power- 
law scalings we label, in Fig.|6] various compensated spectra 
as follows: 

label E: E{k)e^^^h^'\ (6) 



label F : 



(7) 
(8) 
(9) 



\F{k)\e-^^'\"'/2kf, 

label M : M(k)kf(k/k,y^^^ /M, 

label H : \H(k)\kfe'^^'^k^'^E/2M, 

where k^ - J kM(k) dk/M is the wavenumber where the mag- 
netic energy spectrum peaks and ex is the kinetic energy dis- 
sipation per unit mass. The compensated kinetic energy spec- 
trum shows a bottleneck that is clearly stronger than in the 
case without helicity (e.g., Kaneda et al. 2003, Haugen & 
Brandenburg 2006). The kinetic helicity spectrum shows a 
similar spectrum that also has a strong bottleneck, which is 
particularly evident when it is compensated by k^^^. The 
existence of a k^^^ subrange for the modulus of the kinetic 
helicity spectrum is well known from early closure calcula- 
tions (Andre & Lesieur 1977), and has also been seen in di- 
rect numerical simulations (Borue & Orszag 1997, Branden- 
burg & Subramanian 2005a) and in shell model calculations 



Fig. 7. — Relative spectral magnetic helicity, k\H(k)\/2M(k), compensated 
by (kjkf)"!^, for PrM ranging from 0.01 to 1. 

(Ditlevsen & GiuUani 2001). Such a scaling implies that the 
relative spectral kinetic helicity, 

IlKik) = F(k)/2kE(k), (10) 

decreases toward small scales like A:"' and has a maximum at 
k = kf with KKikf) = 0.96. At that scale, the relative magnetic 
helicity, 

'RM(k) = kH(k)/2M(k), (11) 

is -0.15, -0.08, and +0.42 for PrM = 1, 0.1, and 0.01, respec- 
tively. The realizability condition implies that the moduli of 
I^Kik) and "R/nik) are less than unity (Moffatt 1969). The pos- 
itive sign for Pr^ - 0.01 agrees with the idea that the helical 
driving of the flow imprints a helical field of the same sense 
at the same scale. Owing to an inverse cascade of magnetic 
helicity (Pouquet et al. 1976), H{k) is of opposite sign at large 
scales. While this is very clearly established in the nonlin- 
ear regime (BOl) or for small values of Pim, ^ larger range 
of scales attains negative values during the linear stage when 
PrM = 0.1 and 1. 

For the magnetic helicity we also find an approximate 
\H(k)\ ~ k^^^^ spectrum, which is different from the nonlinear 
case when the current helicity, C(k) = k^H(k) shows a k^^^^ 
spectrum (Brandenburg & Subramanian 2005a). Assuming 
that M(k) ~ k^^^, the relative spectral helicity would seem to 
decrease now more rapidly like k\H{k)\/2M(k) ~ fc"'^^^. Fig- 
ure |7] shows that the correspondingly compensated magnetic 
helicity to energy ratio changes now less strongly in the range 
4 < k/ki < 40. 

3.4. Saturation regime 

Eventually the initial exponential growth comes to a halt 
and is followed by a resistively long saturation phase during 
which a large-scale magnetic field develops at wavenumber 
ki, regardless of the value of kj. Owing to the use of periodic 
boundary conditions, this large-scale field tends to be force- 
free and fully helical, and its energy per unit volume is by a 
factor kf/ki -4 larger than the value at kf, which in turn is 
comparable to the kinetic energy per unit volume. For details 
see BOl . Here we only consider the end of this slow saturation 
phase. Compensated kinetic and magnetic energy spectra are 
shown in Fig.[8]for magnetic Prandtl numbers ranging from 1 
to down to 10"^. 

In the final saturated state, and especially for Pr^ = 1, the 
M(k) and E{k) spectra are nearly on top of each other with 
M(k) being slightly larger than E{k) by 20%, which is qualita- 
tively similar to the non-helical case (cf. Haugen et al. 2003). 
There are indications of a somewhat shallower spectrum due 
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Fig. 8. — Kinetic and magnetic energy spectra in the saturated regime for 
PrM = 10^3 with Re = 4400, Pr^ = 10"^ with Re = 2300, Pr^ = 0.1 with 
Re = 1200, and Pr^ = 1 with Re = 450. All spectra are compensated by 
. The ohmic dissipation wavenumber, k^ = (eAf/?/^)'''', is indicated 
by an arrow. The viscous dissipation wavenumbers are 430, 350, 290, and 
180 for PrM = 10"^ 10"-, 0.1, and 1, respectively. 

to a bottleneck effect both for kinetic and magnetic energies 
just before the two enter the viscous and resistive dissipation 
ranges. Also for Pim = 0.1 there is a short range where M{k) 
exceeds E{k), but then, not surprisingly, M{k) turns into the 
dissipation range before E{k) does. This is even more clearly 
the case for Pim = 0.01. 

Low-PrM turbulence has the interesting property that for 
given numerical resolution much larger fluid Reynolds num- 
bers can be achieved than for Pr^ - 1 . This is simply because 
almost all the energy is dissipated resistively, and the energy 
that continues along the kinetic energy cascade is compara- 
tively weak, so not much viscosity is needed for dissipating 
the remaining kinetic energy. In fact, for Pr^ = 0.01 we were 
able to go to Re - 2300 with a resolution of only 512^ mesh- 
points. For PrM = 0.1 and 1 and the same resolution we could 
only go to Re = 1200 and 450, respectively. 

3.5. Diverting most of the energy into Joule heat 
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Fig. 9. — Dependence of the fractional kinetic and magnetic energy dissipa- 
tion rates. Note that the fractional kinetic energy dissipation decreases with 
decreasing Pr^ to the 1/2 power 

As has recently been stressed by Mininni (2007), an in- 
creasing fraction of energy is being dissipated via Joule dis- 
sipation, as PrM decreases. In Fig.|9] we plot the dependence 
of the kinetic and magnetic energy dissipation rates per unit 
mass, eK - {2pvS^)/po and e^f = {WoJ^)/po, relative to the 
total dissipation, = ex + sm, versus Pr^f- The data are well 
described by a power-law fit of the form 



0.37Pr;^l 



(12) 



Thus, for Pr^f = 1 about the 37% of the energy is dissipated 
into viscous heat, while 63% is dissipated via Joule dissipa- 
tion. This is similar to the case of non-helical hydromag- 
netic turbulence (Haugen et al. 2003), where these numbers 
are about 30% and 70%, respectively. 

In turbulence the energy dissipation is generally propor- 
tional to U^/L, where U is the typical velocity and L is a 
typical length scale. Conventionally one defines a dimension- 
less dissipation parameter as 



(13) 



where U is the one-dimensional rms velocity, which is related 
to Mi-ms via = Mnn.s/3i ^nd L is the integral scale and is re- 
lated to kf via jTilkf. In non-helical turbulence this value is 
typically around 0.5 (see also Pearson et al. 2004), but this 
value has never been determined for hydromagnetic turbu- 
lence with helicity. An exception is the work of Blackman & 
Field (2008), who considered a range of power-law scalings 
for kinetic and magnetic energy spectra to calculate analyti- 
cally the dissipation rates. 

It turns out that for our runs, ~ 1.5, i.e. ^ 3 times larger 
than the usual value; see Fig.[TOl Let us now discuss possible 
reasons for this difference. In the definition of the quantity 
one assumes that the energy flux scales with /L. However, 
U is based on the typical rms velocity. In the presence of a 
strong dynamo-generated magnetic field it may be sensible to 
base it on a combination of typical velocity and magnetic field 
strength. In our case we have {B^ I hq) I {pU^) ^ 2, so f/ would 
need to be scaled up by a factor V3, which reduces by a 
factor Ir'l- a! 5 to about 0.3. This value is nearly independent 
of the value of Pr^, which was also found by Blackman & 
Field (2008) under plausible assumptions. 

3.6. Helicity spectra 

As mentioned before, in the nonlinear regime both kinetic 
and current helicities, F(k) and C(k) - k^H{k), respectively. 
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Fig. 10. — Dimensionless total energy dissipation rate, C^, as a function of 
Prj^. Error bars have been estimated based on averages taken over each third 
of the full time series. 
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Fig. 1 1 . — Spectral cun'ent and kinetic helicity ratios for the same four 
runs shown in Fig. [8] The upper panel is for Pr^ = 1 (thin line) and 0.1 
(thick line), while the lower panel is for 10"^ (thin line) and 10"' (thick line). 
Note that for Pr^ = 1 and 0.1 the profiles of {klkf)\<RMiK)\ and (klkf)\'RK(k)\ 
are reasonably flat in the range 6 < kjk\ < 14. The 1/2 slope is shown for 
comparison. 

are expected to display a forward cascade with a k^'^^'^ spec- 
trum. If this is true, we would expect that within some 
wavenumber interval \'Ri!^(k)\ and decrease with in- 

creasing k like ^ '. In Fig.[TT]we show the correspondingly 
compensated relative kinetic and magnetic helicity spectra. 
It turns out that they are surprisingly similar regardless of 
the value of Pr^- For Pr^ = 1 and 0.1 the compensated 
profiles of 'Rui^ and %k{^ are reasonably flat in the range 
6 < kikx < 14. However, for PrM - 10"^ and 10"^ the com- 
pensated profiles show an increase proportional to A:''^. The 
fact that the anticipated fc"' scaling occurs only for magnetic 
Prandtl numbers down to 0. 1 and only over an extremely short 
range may indicate that our Reynolds numbers are still too 
small to yield conclusive results. Especially at smaller scales, 
and certainly in the runs with the smallest Pr^, the compen- 



sated relative kinetic and magnetic helicity spectra are com- 
patible with a fc'^^ slope. This would imply a k"''^^ spectrum 
for the kinetic and current helicities, which is shallower than 
that anticipated for a forward cascade, but still steeper than 
that in the case of equipartition. We emphasize again that this 
applies to the resistively controlled regime. 

In Fig.[TT]we see that at k = kf both JIk and Km are close 
to unity. This indicates that velocity and magnetic fields are 
nearly fully helical. However, for k > kf the velocity and 
magnetic fields become less helical, because the compensated 
relative helicities in Fig. [TT] increase with k not faster than to 
the 1/2 power. On the other hand, for k - k\ the magnetic field 
is again fully helical, i.e. is equal to kijkf - 

1 1 A, but its helicity has the opposite sign. 

In Table[T]we compare the values of Hhiik) during the lin- 
ear and nonlinear stages at the wavenumbers k\ and kf for the 
three or four values of Pr^- Note that during the nonlinear 
stage 'RM(ki) and are of opposite sign. The former 

is close to -1 while the latter increases from 0.52 to 0.72 as 
PrM decreases. As already indicated in ^3.31 during the linear 
stage, the two are of opposite sign only for Pr^ = 0.01, while 
for larger values of Pr^ a larger range of scales appears to 
be affected by the inverse transfer of magnetic helicity caus- 
ing KMikf) to be negative. It would be tempting to try and 
model this behavior using, for example, the four-scale helical 
dynamo model of Blackman (2003). 

TABLE 1 

Comparison of !?^m(^i ) and KMikf) during the linear and nonlinear stages 

FOR different VALUES OF PTjK . 





linear 


nonlinear 


PrM 


ki 


kr 


ki 


kf 


10-3 






-0.993 


+0.72 


10-2 


-0.88 


+0.42 


-0.994 


+0.66 


10-' 


-0.59 


-0.08 


-0.993 


+0.59 


1 


-0.41 


-0.15 


-0.993 


+0.52 



3.7. Effects on the velocity pattern 

In Fig. [12] we compare visualizations of B- and U- for all 
four values of Pr^- The velocity and magnetic field pat- 
terns are surprisingly similar for all four values of Pr^- Only 
for PrM - 10-^ and 10"^ the magnetic field appears notice- 
ably smoother than in the other two cases. The velocity field 
shows a marked anisotropy with small-scale elongated pat- 
terns aligned with the local direction of the mean magnetic 
field, which is here of the form B ~ (0, sinA:iJic, - cos ^ix). 
The anisotropy in the velocity can still be seen for small val- 
ues of PrM, but the small-scale patterns are slightly smoother 

4. CONCLUSIONS 

In many astrophysical bodies the magnetic Prandtl number 
is small, while in most simulations its value is chosen to be 
close to unity. As we have shown here, this mismatch is of 
relatively minor consequence for large-scale dynamos that are 
driven by helical forcing. 

In the nonlinear stage, the velocity and magnetic field pat- 
terns are remarkably independent of the value of PrM- The 
only thing that changes is the length of the inertial range. A 
small magnetic Prandtl number simply means that the mag- 
netic energy spectrum turns into the dissipation range more 
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Fig. 12. — Visualizations of B. and U- for Pr^ = 10 at Re : 
(right). The orientation of the axes is the same as in Fig.[T] 



4400 (left), Pru = 10"^ at Re = 2300, Pr^ = 0.1 at Re = 1200, and Pr^ = 1 at Re = 450 



quickly than the kinetic energy spectrum. It also means that 
essentially all the energy is dissipated via Joule heat. This was 
recently also demonstrated by Mininni (2007). One reason is 
that the case of fully helical turbulence studied in the present 
paper is a particularly simple one, because it leads to uniform 
mean-field dynamo action with large-scale pattern formation 
covering the entire domain. In this paper we have seen that, at 
least for values of Re up to 4400, the dynamics of this large- 
scale pattern, i.e., of the large-scale magnetic field, is quite 
independent of how long the inertial range of the turbulence 
is. In the absence of helicity, there is only small-scale dynamo 
action, which is driven by the dynamics at the smallest pos- 
sible scale, i.e. the resistive scale. In that case it does matter 
what the dynamics of the turbulence is at that scale. How- 
ever, in that case it has not yet been possible to find dynamo 
action for values of Pr^ down to the values considered here. 
Nevertheless, it is possible that even in that case there is an 
asymptotic regime for large enough values of Rcm where the 
dynamics of the magnetic field is independent of the value of 
PrM, even though this regime is not yet accessible with present 
day computers. 

To estimate the value of the magnetic Prandtl number in 
dense astrophysical bodies, one has to use the Spitzer formu- 
lae for T] and v. The resulting magnetic Prandtl number is 
(e.g., Brandenburg & Subramanian 2005b) 

Pr,.l.lxlO-^f^)7 P Xi—V^ (14) 
llO^K/ \0.1gcm-3/ \ 20 / ^ ' 

so at the bottom of the solar convection zone the magnetic 
Prandtl number is clearly rather small (~ 10"**). Neverthe- 
less, simulations of solar and stellar dynamos available so far 
PrM are set to values of the order or unity. Although we have 
shown here that the resulting large-scale fields are similar to 



the more realistic case of small values of Pr^, an important 
difference is that the small-scale dynamo may be more pro- 
nounced when Pr^ is of order unity. In practice this means 
that a positive detection of dynamo action in a simulation 
might not necessarily be relevant for understanding the Sun, 
unless suitable conditions for the excitation of large-scale dy- 
namo action are also met. On the other hand, once the large- 
scale dynamo is really excited, and if it is fully saturated, it is 
then quite feasible to lower the value of Pr^ significantly — 
without losing the large-scale dynamo. In fact, lowering Pr^ 
in a saturated large-scale dynamo means that most of the en- 
ergy will be dissipated via Joule heating, and that the kinetic 
energy cascade only carries a small fraction of the total en- 
ergy. This allows us to increase the value of Re, and hence 
to decrease the viscosity and thereby the value of Pr^ even 
further. Simulations of large-scale dynamo action in turbulent 
convection (Kapyla et al. 2008) provide one example where it 
is indeed feasible to lower Pr^, although in that case the sys- 
tem is not uniform and so energy dissipation via Joule heating 
is only possible in those locations where the dynamo is strong 
enough (P. J. Kapyla 2008, private communication). 
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